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We report a measurement of the i/^^-nucleus inclusive charged current cross section on iron 

using data from exposed to the J-PARC neutrino beam. The detector consists of 14 modules in 
total, which are spread over a range of off-axis angles from 0° to 1.1°. The variation in the neutrino 
energy spectrum as a function of the off-axis angle, combined with event topology information, is 
used to calculate this cross section as a function of neutrino energy. The cross section is measured to 
becr°°(l.l GeV) = 1.10±0.15 (10“^®cm^/nucleon), cr°°(2.0 GeV) = 2.07±0.27 (10/nucleon), 
and cr°°(3.3 GeV) = 2.29 ± 0.45 (10/nucleon), at energies of 1.1, 2.0, and 3.3 GeV, respec¬ 
tively. These results are consistent with the cross section calculated by the neutrino interaction 



generators currently used by T2K. More importantly, the method described here opens up a new 
way to determine the energy dependence of neutrino-nucleus cross sections. 
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I. INTRODUCTION 

Many recent long baseline neutrino oscillation exper¬ 
iments use muon neutrino beams, with neutrino ener¬ 
gies ranging from sub-GeV to a few GeV. The observed 
neutrino-nucleus charged-current (GG) interactions are 
then used to infer neutrino oscillation parameters. In this 
energy region, GG quasi-elastic and GG single-pion pro¬ 
duction reactions dominate the total cross section, and 
so understanding these channels is essential for precision 
measurements of the oscillation parameters. Measure¬ 
ments of exclusive cross sections, however, are compli¬ 
cated by re-interactions of the final state hadrons as they 
exit the nucleus, known as final state interactions (FSI). 
FSI can absorb or produce particles, resulting in a dif¬ 
ferent set of particles entering the detector than would 
be expected from the initial interaction. For example, 
the pion from a GG single-pion interaction might be ab¬ 
sorbed in the nucleus, so the observable final state is 
similar to that from a GG quasi-elastic event. The GG 
inclusive channel is much less sensitive to these effects, 
since it only requires the detection of a charged lepton 
(muon) from the interaction. A precise measurement of 
this channel, combined with exclusive measurements, will 
help improve our understanding of neutrino interactions 
in this energy region. 

So far, the MINOS and T2K experiments have mea¬ 
sured the GG inclusive cross section on iron [1, 2] us¬ 
ing neutrino beams which cover the few-GeV region. The 
former experiment measured the GG inclusive cross sec¬ 
tion for neutrinos with energies above 3.5 GeV using the 
MINOS near detector. The latter used the T2K near de¬ 
tector, INGRID, to measure a flux-averaged GG inclusive 
cross section at a mean energy of 1.51 GeV, where the 
r.m.s spread of the neutrino energy was 0.76 GeV (0.84 
GeV) below (above) this mean energy. The MINERz/A 
experiment also measured the GG inclusive cross section 
on iron, but only the cross-section ratio of iron to GH 
has been published [3]. The GG inclusive cross section 
on iron has not yet been measured in the 2-3 GeV energy 
range, and a measurement covering 1-3 GeV would pro¬ 
vide a consistency check between the T2K and MINOS 
results. 

The T2K experiment is a long baseline neutrino oscil¬ 
lation experiment in Japan [4]. T2K utilizes an almost 
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pure beam produced as the decay product of tt+’s and 
iC+’s. The beam is first measured by the near detector, 
ND280, located 280 m downstream from the pion pro¬ 
duction target. After traveling 295 km, the neutrinos are 
then observed at the far detector, Super-Kamiokande. 
Oscillation parameters are determined by comparing the 
neutrino interactions observed at the near and far detec¬ 
tors. 

T2K was built around the “off-axis beam method”, 
where detectors are intentionally placed off the central 
axis of the neutrino beam (hereafter beam-axis). The 
angle with respect to the beam-axis is called the off-axis 
angle: ^OA- The direction of the neutrino parent particles 
is distributed around the beam-axis, so 6 >oa is approxi¬ 
mately equal to the angle between the parent particle and 
neutrino directions. In this case, the energy of a neutrino 
produced from the two-body decay (tt ^ + ji) can be 

expressed as follows: 

9 9 

—m,, 

E. = ’ ( 1 ) 

2(-£/7r Ptt cos 0q\) 

where and are the masses of the pion and muon 
whilst and are the energy and momentum of the 
pion. The relation between Ej^ and for different ^oas 
is plotted in Fig. 1, showing the maximum neutrino en¬ 
ergy reducing as Oqa increases. This indicates that the 
energy spectrum of the neutrino beam peaks at a lower 
energy and has a narrower width as Oqa increases. In 
the T2K experiment, Oqa = 2.5° was chosen so that the 
neutrino flux peaks around 0.6 GeV, an energy which 
maximizes the oscillation probability of the muon neu¬ 
trino at the far detector. 



(GeV/c) 


FIG. 1. Neutrino energy as a function of the pion momentum 
for different Oqas. 

The T2K INGRID detector is installed on the beam- 
axis at the near site. It consists of 14 identical modules, 
which are spread over a range of Oqa from 0° to 1.1°. 
Thus, the peak of the neutrino energy spectrum differs 
among the modules as in Eq. (1). 
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In this paper we present a measurement of the in¬ 
clusive CC cross section on iron in the energy range of 
1-3 GeV with INGRID. This analysis uses data collected 
from 2010 to 2013, corresponding to 6.27 x 10^^ protons 
on target (p.o.t.). The neutrino interactions at different 
INGRID modules, which are distributed at different posi¬ 
tions and thus observed different beam spectra, is used to 
extract the energy dependence of the cross section. The 
topology of each event, which is based on the kinematics 
of the outgoing muon, is also used to further improve the 
sensitivity of this measurement to the neutrino energy, 
since the two are directly related. The different energy 
spectra and event topologies are combined to construct 
a probability density function (PDF), which is used to 
measure the cross section using the least method. 

The paper is organized as follows. In Sec. II, we de¬ 
scribe the T2K near detector, INGRID. Section III intro¬ 
duces the Monte Garlo (MG) simulation used to predict 
neutrino event rates at the INGRID detector and de¬ 
scribes the systematic uncertainties associated with this 
model. Next, the analysis method used to extract the 
energy dependence of the cross section is explained in 
Sec. IV with a discussion of the remaining systematics in 
Sec. V. Finally, the result of the analysis is presented in 
Sec. VI. 


II. THE T2K NEAR DETECTOR: INGRID 

Situated 280 m downstream from the pion production 
target, the INGRID detector monitors the neutrino beam 
direction and intensity. It consists of 14 identical mod¬ 
ules, each of which is composed of 9 iron target plates and 
11 scintillator tracking planes. The iron plates and the 
tracking planes are stacked in alternating layers forming 
a sandwich structure, as shown in Fig. 2. 


Veto planes 



FIG. 2. Schematic drawing of an INGRID module. Each 
module has 9 iron target plates and 11 tracking planes, with 
4 veto planes covering the side surfaces. 


from outside of the module. Each tracking plane has 
two layers of scintillator bars aligned orthogonally to one 
another, enabling particles to be tracked in all 3 dimen¬ 
sions as they pass through the plane. The veto planes are 
also formed from scintillator bars. The bars are coated in 
Ti 02 reflectors to help contain scintillation light, which is 
then captured by wavelength shifting (WLS) fibers which 
run through the center of the bars. This light is then read 
out by a Multi-Pixel Photon Gounter (MPPC) [5, 6] and 
the resultant signal is digitized by the Trip-t front-end 
board (TFB) [7], to give the integrated charge and tim¬ 
ing information. The integration cycle of the electronics 
is synchronized with the neutrino beam pulse structure, 
ensuring all data is captured. 

The modules are installed in a cross shape centered on 
the beam-axis. An overview of the INGRID detector is 
shown on the top in Fig. 3. An ID is assigned to each 
module as shown on the bottom in Fig 3. Further details 
of the detector and the basic performance of INGRID can 
be found in Ref. [8]. 



Each of the iron plates has dimensions of 124 x 124 cm^ 
and a thickness of 6.5 cm, providing a total iron mass of 
7.1 ton per module. The module is surrounded by scintil¬ 
lator veto planes, which detect charged particles coming 


EIG. 3. Top: The INGRID detector. The modules are ar¬ 
ranged inalOmxlOm cross. The two off-axis modules not 
located in the arms of the cross are not used in this analysis. 
Bottom: Module ID given to each module. 
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III. SIMULATING NEUTRINO INTERACTIONS 
IN THE INGRID DETECTOR 

In this analysis the cross section is measured by com¬ 
paring data to a prediction, which is calculated using 
three sequential MC simulations: 

1. Prediction of the neutrino flux at the INGRID mod¬ 
ules. 

2. Simulation of the neutrino-nucleus interactions in 
the iron target (^^Fe). 

3. Propagation of final state particles through the de¬ 
tector and modeling of its response. 

After these processes, an event selection, which is de¬ 
tailed in Sec. IV B, is applied to the output in the same 
way as the data. 

A. Neutrino flux 

1. Flux prediction 

A detailed description of the neutrino flux predic¬ 
tion can be found in Ref. [9]. In the simulation, pro¬ 
tons are impinged upon the carbon target to produce 
hadrons, which decay into neutrinos. FLUKA2008 [10] 
and GEANT/GCALOR [11] are used to model hadron 
production in the target and surrounding material. Prop¬ 
agation of the resulting particles through the electromag¬ 
netic horns, which focus the charged hadrons along the 
beam-axis, is simulated using dedicated, GEANT3 [12]- 
based code, which also models the subsequent decay of 
the particles. Eor each hadron decay mode which pro¬ 
duces neutrinos, the probability of the neutrinos to be 
emitted in the direction of the INGRID detector is cal¬ 
culated. The flux prediction is obtained by weighting the 
generated neutrinos with these probabilities. The flux is 
then tuned using hadron interaction data, primarily from 
the NA61/SHINE experiment [13]. Other hadron pro¬ 
duction data (Eichten et al [14] and Allaby et al [15]) 
are also used to tune the simulation in regions of the 
hadron production phase space that are not covered by 
the current NA61/SHINE measurement. In this analysis, 
the NA61/SHINE data taken in 2007 is used to correct 
the neutrino flux [16, 17]. 

The neutrino flavor content across different energy re¬ 
gions at module 3 (one of the center modules) is sum¬ 
marized in Table I. This shows that muon neutrinos ac¬ 
count for >^95% of the total flux for Ejy < 3 GeV. The 
muon neutrino flux fraction then falls to less than 90% 
for El, >3 GeV, where + z^e account for >^10% of the 
total flux. The flavor content of the neutrino flux at the 
other modules is similar to that at module 3. 

Eigure 4 shows the obtained muon neutrino flux spec¬ 
tra at the INGRID modules. The neutrino energy spec¬ 
trum changes with module position, with the spectrum 


TABLE I. Fraction of the integrated flux by neutrino flavor 
in each energy range at module 3. 



Neutrino energy range (G 

eV) 

Flavor 

0-1 

1-2 

2-3 

3-4 

>4 


94.2% 

96.8% 

95.4% 

89.7% 

86.5% 


4.8% 

2.7% 

3.8% 

7.9% 

9.3% 


0.9% 

0.5% 

0.7% 

2.0% 

3.5% 

Ve 

0.1% 

0.0% 

0.1% 

0.3% 

0.6% 


at module 0 softer than that at module 3. This is be¬ 
cause module 0 is located at Oqa =1.2° from the neu¬ 
trino beam-axis, and so the neutrino flux passing through 
it is shifted to lower energies due to the off-axis beam 
effect. This feature is, indeed, essential for this cross- 
section measurement. 


xlO'^ 



FIG. 4. Muon neutrino fluxes at modules 0, 1, 2 and 3. 


2. Flux uncertainties 

The systematic error on the neutrino flux prediction 
comes from uncertainties in hadron production and from 
errors in the measurement of the proton beam, the horn 
current, and the target alignment. The uncertainties 
in the hadron production are mainly driven by uncer¬ 
tainties in the NA61/SHINE measurements.and those in 
measurements by Eichten et al. and Allaby et al. The 
second category of flux errors is associated with inher¬ 
ent uncertainties and operational variations in the beam¬ 
line conditions. They include uncertainties in the proton 
beam position, the beam direction, the absolute horn cur¬ 
rent, the horn angular alignment, the horn field asymme¬ 
try, the target alignment and the proton beam intensity. 
The method used to estimate these flux uncertainties is 
described in Ref. [9]. Eigure 5 shows the flux error at 
module 3, which includes all sources of uncertainty, and 
demonstrates that the systematic error on the neutrino 
flux is dominated by the uncertainties in the hadron pro¬ 
duction model. The propagation of the flux error in this 
analysis is described in Sec. V A. 
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FIG. 5. Fractional flux error at module 3, including all 
sources of uncertainty. 


B. Neutrino-nucleus interaction simulation 

Neutrino-nucleus (Fe) interactions in the INGRID 
module are simulated by a neutrino event generator, 
which is a composite of different neutrino interaction 
models. NEUT (ver.5.1.4.2) [18] is used as the primary 
event generator in this analysis. GENIE (ver.2.8.0) [19], 
a different neutrino interaction simulation package, is 
also used for comparison. This section first describes the 
various interaction models simulated in the NEUT event 
generator, then explains the systematic uncertainties as¬ 
sociated with each model. Details of the event generators 
used in T2K can be found in Ref. [20] . 


1. The NEUT neutrino event generator 

Given a neutrino energy and a detector geometry, 
NEUT determines the interaction mode of an event and 
calculates the kinematics of the final state particles. It 
also simulates ESI for hadrons as they traverse the target 
nucleus. The following interaction modes are provided 
for both charged-current (CC) and neutral-current (NC) 
interactions by NEUT: 

• Quasi-elastic scattering (CCQE or NCQE) 

• Resonant pion production (CCItt or NCItt) 

• Deep inelastic scattering (CCDIS or NCDIS) 

• Coherent pion production 

Here N and N' denote nucleons, I is the lepton, and A 
is the target nucleus. 


Neutrino interaetion model uncertainties 

Table II summarizes the parameters used for modeling 
neutrino interactions in NEUT. The systematic param¬ 
eters listed in the table were evaluated in the previous 
analyses of neutrino oscillation from T2K [21, 22] and 
fall into the following four categories: 


, and the nuclear model: 

These parameters are used for modeling CCQE and 
CCItt interactions. A 20% error is assigned to 
the axial vector masses, which comes from a com¬ 
parison of external measurements of these parame¬ 
ters. The uncertainty on the Eermi momentum and 
binding energy are estimated from electron scat¬ 
tering data [23]. The uncertainty in the nuclear 
model is evaluated by exchanging the REG nuclear 
model with the spectral function model described 
in Ref. [24]. 

TT-less A decay and W shape: 

In the resonant pion production process, baryon 
resonances, mainly As, can interact with other nu¬ 
cleons in the target nucleus, and disappear without 
pion emission. The tt-Icss A decay parameter is in¬ 
troduced to take into account the uncertainties in 
this process. The W shape parameter is introduced 
to modify the shape of the momentum distribution 
of pions from NC single pion production interac¬ 
tions so that it matches MinibooNE data [25]. 

Normalization parameters: 

Normalization parameters are used to change the 
overall normalization of the cross section. The nor¬ 
malizations for CCQE and CCItt are defined sepa¬ 
rately for different energy regions. The uncertain¬ 
ties on the normalizations are mostly determined 
from the MiniBooNE data. The CC other shape 
parameter is introduced as an energy dependent 
uncertainty on CC DIS and CC resonant interac¬ 
tions, where the resonance decays into a nucleon 
and photon, kaon or eta. According to the MINOS 
measurement [1], the relative uncertainty on the 
CC inclusive cross section on iron, which is domi¬ 
nated by CCDIS, is approximately 10% at 4 GeV. 
Using this as a reference point, the error on the 
CCDIS and CC resonant cross section is scaled us¬ 
ing the following formula: 

^^CCother _ 0-4 GcV , . 

^CCother (GcV) 

Although the error goes to infinity as Ej, ap¬ 
proaches 0 GeV, this error does not have a signifi¬ 
cant contribution to the total cross section error for 
the lower neutrino energy region because the inter¬ 
actions have a threshold energy of approximately 
0.6 GeV and a small cross section in the 1 GeV 
energy region. Einally, the Itt Ej^ shape parameter 
is a weighting factor as a function of neutrino en¬ 
ergy. This is introduced to cover the discrepancy 
between the MiniBooNE measurement of the CCItt 
cross section versus E^^ and the NEUT prediction 
using the best fit parameters obtained from the fit 
to the MinibooNE data [26]. This discrepancy is 
as large as 50% at 600 MeV, so an error of 50% is 
assigned to this weighting factor. In the nominal 
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TABLE IL Neutrino interaction systematic parameters, nominal values, uncertainties (Icr), and interaction types (CC, NC 
or CC+NC). The first, second, and third groups represent the model parameters, the ad hoc parameters, and the pion ESI 
parameters respectively. A 1 or 0 in the nominal value column means that the effect of the systematic parameter is implemented 
or not implemented by default [21, 22]. 


Parameter 

Nominal value 

Uncertainty (Icr) 

Interaction type 

mT 

1.21 GeV 

0.20 GeV 

GG 

MREs 

1.21 GeV 

0.20 GeV 

GG+NG 

Fermi momentum (Fe) 

250 MeV/c 

30 MeV/c 

GG 

Binding energy (Fe) 

33 MeV 

9 MeV/c 

GG 

Spectral function 

0 (off) 

1 (on) 

GG 

TT-less A decay 

0.2 

0.2 

GG + NG 

W shape 

87.7 MeV 

45.3 MeV 

GG+NG 

CCQE normalization <1.5 GeV) 

1 

0.11 

GG 

CCQE normalization (1.5 < Eu < 3.5 GeV) 

1 

0.30 

GG 

GGQE normalization (Eu > 3.5 GeV) 

1 

0.30 

GG 

GGItt normalization (Eu < 2.5 GeV) 

1 

0.21 

GG 

GGItt normalization (Eu > 2.5 GeV) 

1 

0.21 

GG 

GG coherent normalization 

1 

1.0 

GG 

GG other shape 

0 

0.1 at P^=4.0 GeV 

GG 

NG Itt^ normalization 

1 

0.31 

NG 

NG coherent pi normalization 

1 

0.30 

NG 

NG Itt^ normalization 

1 

0.30 

NG 

NG other normalization 

1 

0.30 

NG 

Itt Eu shape 

0 (off) 

0.50 

GG + NG 

Pion absorption 

1 

0.5 

GG+NG 

Pion charge exchange (Pj, < 500 MeV/c) 

1 

0.5 

GG+NG 

Pion charge exchange (P^r > 400 MeV/c) 

1 

0.3 

GG+NG 

Pion QE scattering (P^r < 500 MeV/c) 

1 

0.5 

GG+NG 

Pion QE scattering (Ptt > 400 MeV/c) 

1 

0.3 

GG+NG 

Pion inelastic scattering 

1 

0.5 

GG+NG 


NEUT, this weighting is not applied (“off” as in 
Table II). 

Pion FSI: 

There are three pion FSI processes of interest in the 
T2K energy range: absorption, charge exchange 
and QE scattering. In addition to these interac¬ 
tions, the particle production process, defined as 
“inelastic scattering”, was considered, since it is the 
dominant process at higher pion energies. Uncer¬ 
tainties on the FSI parameters are estimated using 
external data sets [27]. 

Propagation of these uncertainties is described in 
Sec. VB. 


C. Detector simulation 

The particle type and kinematic information provided 
by NEUT is passed to the detector simulation built 
within a GEANT4 framework [28] . All the detector com¬ 
ponents are modeled in the simulation. The energy de¬ 
posited by each particle in the scintillator planes is con¬ 
verted into a number of photo-electrons, taking into ac¬ 
count the non-linear response of the scintillator, light col¬ 
lection efficiency and attenuation of the WLS fiber, and 
the non-linearity of the MPPC response. The non-linear 


response of the ADCs on the front-end electronics is also 
modeled based on the results of charge injection tests. 

Particles generated in the wall upstream of the IN¬ 
GRID detectors are also propagated into the detector 
simulation, and are treated as a background (BG) source. 

The hadronic interaction of particles in the detector is 
simulated by GEANT4 using the FTFP_BERT physics 
list. In this physics list, the GEANT4 Bertini cascade 
model is used to simulate nuclear reactions by hadrons 
with kinetic energies below 5.5 GeV. For particles with 
kinetic energies above 5 GeV, the list uses the Fritiof 
model [29, 30]. 


IV. ANALYSIS METHOD 
A. Overview 

First, neutrino interactions are selected in each IN¬ 
GRID module. The different neutrino energy spectra 
sampled by the different modules provide a way of ex¬ 
tracting the energy dependent cross section. The selected 
events are further categorized according to the topology 
of their final-state muon in order to improve the sensitiv¬ 
ity of the samples to the energy of the incoming neutrino. 
A PDF is then constructed relating the different INGRID 
modules and event topologies to the neutrino energy. Fi¬ 
nally, to extract the GG inclusive cross section, a fit 









is performed between the selected events and the PDFs. 


B. Neutrino event selection 

A detailed description of the event selection for neu¬ 
trino interactions at INGRID can be found in Refs. [2, 
31]. A brief explanation of each step of the selection is 
given here. A typical selected muon-neutrino interaction 
candidate is shown in Fig. 6. 


5. Timing cut: 

The T2K neutrino beam has eight bunches in 
each beam spill, and each bunch has a width of 
58 nsec. The selected events are required to lie 
within 100 nsec of the expected time of each bunch. 

6. Veto cut: 

The reconstructed track is extrapolated to up¬ 
stream positions in the side veto planes and up¬ 
stream veto planes, and the event is cut if any hits 
are found near the extrapolated entry point. 


Side view Top view 



FIG. 6. Event display of a muon-neutrino event candidate. 
Circles and solid lines represent the hits at the scintillator 
planes and reconstructed tracks, respectively. The neutrino 
beam enters from the left side. 


1. Pre-selection: 

The integrated charge and timing of hits in each 
channel are recorded with a 2.5 photo-electron 
threshold. If there are more than 3 hit channels 
within a 100 nsec time window, these hits are com¬ 
bined to form a hit cluster. The scintillator planes 
are then searched to find those with at least one 
hit in both their X and Y oriented layers. Such a 
plane is called an “active plane”, and events with 
2 or more active planes are selected. 

2. 2D track reconstruction: 

A “cellular automaton” tracking algorithm [32] is 
applied to hits in the X and Y planes to obtain 
tracks in the XZ view and YZ view respectively. 

3. 3D track matching: 

The difference in the most upstream layer hit be¬ 
tween any two tracks in the XZ view and YZ view is 
used to determine if they originate from the same 
vertex. If this difference is greater than 2 planes 
then the tracks are not matched. 

4. Vertexing: 

The vertex is defined as the most upstream hit of 
the track. If there are two or more 3D tracks, a 
check is performed to see if they originate from the 
same vertex or not. 


7. Fiducial volume cut: 

The fiducial volume (FV) is defined as a cube with 
a (±50) x(±50) cm^ transverse area, correspond¬ 
ing to the scintillator bars from the 3rd to 22nd 
channel in the X and Y direction, and from track¬ 
ing plancT^l to 8 in the Z direction (see Fig. 7). 
Events with a vertex inside the FV are selected. 


Tracking plane (X,Y) 





FV 









10cm 100cm 10cm 


FIG. 7. Fiducial volume of the INGRID module viewed from 
of the side (top) and front (bottom). 


After the event selection above, corrections are applied 
to the MG to account for differences in the individual iron 
target masses, the observed background rate, and the 
number of dead channels. We also apply a correction for 
event pileup, which depends on the beam intensity and 
results in a loss of efficiency [31]. Table III summarizes 
the number of events predicted by the MG simulation 
and observed in the data. 

The selected events contain the muon-neutrino signal 
as well as background events, such as and inter- 
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TABLE III. Summary of the neutrino event selection. The p.o.t. used in the MC simulation is normalized to correspond to 
the data set used in this analysis, i.e. 6.27 x 10^°. 



Data 



MC 

beam-related BG 

Total 

Vertexing 
Timing cut 

Veto cut 

FV cut 

3.993 X 10^ 
3.992 X 10’’ 
1.725 X 10’’ 
1.103 X 10^ 

1.655 X 10’’ 
1.655 X 10^ 
1.458 X 10^ 
1.098 X lO’’ 

0.039 X 10’’ 
0.039 X 10’’ 
0.036 X 10’’ 
0.027 X lO’^ 

2.294 X lO’’ 
2.294 X lO’’ 
0.239 X lO’’ 
0.006 X lO’’ 

3.987 X lO’’ 
3.987 X lO’’ 
1.733 X lO’’ 
1.131 X lO’’ 


actions. The other backgrounds come from muons, neu¬ 
trons and photons generated by neutrino interactions out 
of the FV or in the pit wall upstream of INGRID (here¬ 
after called beam-related BG). Since the contamination 
of Pe is negligible, it is not counted in the MC. The angu¬ 
lar distribution of the lepton track for the selected events 
is shown in Fig. 8. Since the vertex is defined as the most 
upstream hit of the track, the angular acceptance is lim¬ 
ited to between 0° and 90°. 



Track angle (degree) 


FIG. 8. Distribution of the reconstructed track angle with 
respect to the beam direction after the event selection. The 
number of events shown in the figure is the total integrated 
over all modules. 

The final selected event sample has >70% efficiency for 
CC interactions from neutrinos with energies >1 GeV, as 
shown in Fig. 9. Around 5% of the selection inefficiency 
in the higher energy region is due to events where the 
muon is produced at a large angle to the Z axis of the 
detector. This results in it escaping the module before 
passing through two active planes. The predicted energy 
spectrum of the reconstructed events at the INGRID 
modules is shown in Fig. 10. 



FIG. 9. Detection efficiency for CC+NC (solid line) and CC 
(dashed line) events. These efficiency curves are estimated 
from the number of events integrated over all modules. 


xlO^ 



FIG. 10. Predicted energy spectrum of the reconstructed 
events at the different INGRID modules. 


1. Downstream (DS-) escaping 


C. Event topology 

To improve the sensitivity of this analysis to the energy 
of the neutrino, the selected events are categorized into 
the following two topologies: 


2. Non-downstream (NonDS-) escaping 

If one of the tracks from the neutrino interaction pene¬ 
trates the most downstream plane, as shown on the left 
in Fig. 11, that event is categorized as DS-escaping. All 
other events, i.e. both side escaping and fully contained 
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events (see the right plot in Fig. 11) are categorized as 
NonDS-escaping. 



> Z 



FIG. 11. Different event topologies. If a track penetrates 
the most downstream plane, the event is categorized as DS- 
escaping (left). The other events are categorized as NonDS- 
escaping (right). 


Events are then further categorized according to the re¬ 
constructed vertex-Z position. Vertex-Z is defined as the 
most upstream active plane number and ranges from 1 
(most upstream) to 8 (most downstream). Events whose 
vertex-Z is in the most downstream plane are greatly af¬ 
fected by uncertainties in the GEANT4 hadron produc¬ 
tion model, so only events with vertex-Z in the range 1-7 
are used in this analysis. 

In total, there are 14 event topologies: 



FIG. 12. Neutrino energy spectra for “NonDS-escaping & 
vertex-Z=7” (solid line) events and “DS-escaping & vertex- 
Z=l” (dashed line) events. The energy spectra for module 0 
are shown, normalized by area. 


• DS-escaping: vertex-Z=l-7 

• NonDS-escaping: vertex-Z=l-7 

Figure 12 shows the energy spectra of “DS-escaping 
& vertex-Z=l” events and “NonDS-escaping & vertex- 
Z=7” events for module 0. The former has a more ener¬ 
getic /i track and are generally produced by higher energy 
neutrinos. The latter, on the other hand, tend to have 
muons produced at a larger angle to the neutrino beam 
or with a lower energy, and so the majority come from 
lower energy neutrinos. 


TABLE IV. Definition of the module groups 


Module Module 

Distance from 

Off-axis Horizontal or 

group 

ID 

the beam axis (cm) 

angle 

Vertical 

1 

0, 6 

450 

1.1° 

Horizontal 

2 

7, 13 

450 

1.1° 

Vertical 

3 

1, 5 

300 

o 

o 

Horizontal 

4 

8, 12 

300 

O 

o 

Vertical 

5 

2,4 

150 

O 

o 

Horizontal 

6 

9, 11 

150 

O 

o 

Vertical 

7 

3, 10 

0 

0° 

(Genter) 


D. Module grouping 

A shift of the neutrino beam direction changes the peak 
of the neutrino energy spectra at the INGRID modules. 
In order to reduce this effect, for the horizontal and ver¬ 
tical directions separately, the two modules at beam-axis 
symmetric positions are grouped together. This results 
in 7 module groups in total, defined in Table IV. 

The number of selected events for each module group 
and each topology is then defined as: 

(3) 

where the indices j and g denote the topology and 
the module group {g = 1,2,..,7), respectively. The 
m and m' indices stand for the module numbers corre¬ 
sponding to each module group. 


E. Detector response uncertainties 

This section introduces two different kind of detector 
response uncertainty: those producing correlation among 
the event topologies and those that do not. They are 
summarized as follow: 

• Uncorrelated error sources 

— Mass of iron plate 
— Pileup correction 

• Correlated error sources 

— Event selection and reconstruction 
— MPPC noise rate 

The treatment of these systematic uncertainties in this 
analysis is described in Sec. V 
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1. Uncorrelated errors 

Iron mass: 

The error on the measurement of the mass of each 
iron plate and the machining tolerance for the plate 
area are taken into account in the systematic error 
on the iron mass. Since these errors are indepen¬ 
dent for each iron plate, an uncorrelated error of 
0.09% assigned to the number of selected events. 

Pileup correction: 

The number of selected events is corrected to ac¬ 
count for event pileup. The correction factor is es¬ 
timated using data sets at different beam inten¬ 
sities. The uncertainty on the correction factor 
comes from the statistical error on the number of 
events, so the uncertainty is uncorrelated between 
the event topologies. An error of 0.5-2.0%, depend¬ 
ing on the event topology, is assigned in this anal¬ 
ysis. 

2. Correlated errors 

MPPC noise: 

MPPC noise hits sometimes result in mis- 
reconstruction of the event vertex or a miscounting 
of the number of active planes, which produces a 
variation in the neutrino event selection efficiency. 
The systematic error caused by the variation in the 
measured noise rate over time is evaluated by alter¬ 
ing the noise rate in the MC. As a result, a 0.1-1% 
error is assigned to the number of selected event in 
each topology. 

Event selection: 

In this analysis, uncertainties in the following event 
selection steps are taken into account: 

• 3D track matching 

• Vertexing 

• Veto cut 

• FV cut 

The systematic error on the number of selected 
events is evaluated by varying the selection thresh¬ 
old for each step, picking the loosest or tightest 
threshold. The change from the nominal threshold 
to the loosest (tightest) is defined as the +lcr(—la) 
change. The resultant fractional variation in the 
number of selected events (= AN/N) due to the 
±la change is computed for both the data and MC. 
Any difference in AN/N between the data and MC 
is then taken as a systematic error. 

Uncertainties on the hit efficiency of the tracking 
planes, the contamination due to beam-related BGs, the 
hit inefficiency of the upstream veto plane, and the track¬ 
ing efficiency were found to be negligible for this analysis 
and are not included in the final result. 


F. Cross-section extraction 

This analysis uses the least method to fit to the 
observed number of events at each module group 
bin: 1-7) and for each event topology bin: 1-14): 

(Xfobs _ tj^cc xfnc X N^^)V 

o ^ J9~ J9~ jgti 

- 

+ J2Mh\VkUAf^ , (4) 

k 

where is the observed number of events, 
and N^^ are the expected numbers of CC events, NC, 
and BG events, and Afk and 14 are the systematic pa¬ 
rameter and the covariance for the error source, re¬ 
spectively. For the covariance term uncertainties in the 
neutrino flux, neutrino interaction model and the detec¬ 
tor response are taken into account and are described 
in Sec. V. The denominator in the statistical term is 
composed of the statistical error on the observed number 
of events the MG statistical error (aTvj^c), and 

the error on the detector response, which is uncorrelated 
among event topologies: (a^det): 



The expected number of GG events in the module 
group and for the event topology is expressed as: 

Xa E (1 + + AX + Afi) 

i 

x0i,-<Tr-e--T , (6) 

• A/^ : systematic parameter for the detector response, 

• A/^^: systematic parameter for the GG interaction model, 

• A/^ : systematic parameter for the flux, 

• (j) \ flux, 

• GG cross section, 

• : detection efficiency for the GG interaction, 

• T : the number of nucleons 

in the fiducial volume of the INGRID module. 

Here i goes over the neutrino energy bins, described in 
Sec. IV G, and Afi is the parameter being fit, which is 
used to represent fractional deviations of the GG inclusive 
cross section. The systematic parameters, A/j^, Afj^ 
and Af!^g^ are also fit to include the effect of these system- 
atics into the cross-section result. The INGRID modules 
are formed from both iron and plastic scintillator (CH). 

The effect on this result coming from the different target 
nucleus for GH interactions is found to be small, so the 
event rate per unit weight on GH is assumed to be equal 
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to that on iron. Afj and Afj^ are systematic parame¬ 
ters representing the uncertainty on the detector response 
and the CC cross-section model respectively for the 
topology bin. These uncertainties change the detection 
efficiency as a function of neutrino energy, resulting in a 
variation in the number of selected events. The difference 
in these uncertainties between the module groups is very 
small, therefore the same parameters are applied to all 
module groups. Affg parameterizes the flux uncertainty, 
changing the normalization of the neutrino flux, in the 
j^th QYieigy bin of the module. The A/j^, Afj^ and 
A fig parameters describe fractional deviations from the 
nominal MC and change the number of events in each 
event topology, module group, and energy bin. 

Since the fraction of NC events in the selected sample 
is very small, the NC events are summed over the entire 
energy region and an averaged flux systematic parameter 
is applied to them: 

We express the number of the NC events as follows: 

E(i+(8) 

i 

The other BG events (the beam related BG and and 
z/g beam flux components) are summed for each module 
and for each topology (Njg). The number of BG events 
in the sample, and their associated errors, are both small, 
so the errors on these BGs are neglected in this analysis. 


G. Energy binning 


The cross section is required to be continuous at the bin 
boundaries and is linearly interpolated between bound¬ 
aries, so Eq. (6) is modified as follows: 


^ s =(1 + A //)-^(1 + A 4 ) 

Li 

xE 


1=0 


l + Afi+ A/»+i A/j _ ^ 


X • T 


(9) 


To interpolate, each energy bin is divided into fine 
bins (/ = 0,1,.., I/^). The energy bin is defined as the 
“global bin” and the energy bin as the “local bin”, 
respectively. The energy range of the global bins and the 
binning for the local bins are summarized in Table V. 
The AfiS are set at 0.5, 0.8, 1.4, 2.6, and 4.0 GeV. 

Finally, the energy dependent GG inclusive cross sec¬ 
tions are extracted as follows. After deriving the fit pa¬ 
rameters using the least method, the cross sections are 
obtained by multiplying those in the original model by 


TABLE V. Summary of the energy range of each global bin, 
the size of each local bin, and the number of the local bins 
per global bin. 

Energy range of Size of each Number of 
global bin (GeV) local bin (MeV) local bins (Li) 


0-0.5 

500 

1 

0.5-0.8 

100 

3 

0.8-1.4 

100 

6 

1.4-2.6 

100 

12 

2.6-4.0 

100 

14 

4.0-30.0 

26000 

1 


^ ^ a/i+a/ 2 ^ ^ ^ A/ 2 +A /3 ^ ^ Ah+AU ^ Taking the 

average of neighboring parameters produces a measure¬ 
ment at the central energy between the bin boundaries, 
since, as in Eq. (9), a linear interpolation is applied be¬ 
tween the neighboring Afi parameters. As a result of this 
averaging, the cross section is measured at 1.1, 2.0, and 
3.3 GeV. The final error on the cross section is smaller 
than those on the Afi parameters due to the negative 
correlations between the Affs. This feature is a result 
of the cross section continuity requirement. 

As seen in Fig. 9, the INGRID detection efficiency for 
GG interactions falls rapidly for neutrinos with energies 
less than 0.5 GeV. Since the event samples are not sen¬ 
sitive to the cross section in this region A/o is not used 
in the final result. For Ej^ > 4.0 GeV there is only a 
small difference in the neutrino energy spectra between 
the INGRID modules. Therefore, sensitivity to the cross 
section for Ey > 4.0 GeV is expected to be worse com¬ 
pared to the lower energy regions. For these reasons it 
was decided, before fitting the data, to use A/1-A/4 to 
measure the cross section at 1.1, 2.0, and 3.3 GeV. 

The PDF of GG events in the global energy binning is 
shown in Fig. 13. Here, the “Fraction” described by the 
z-axis of the figure is obtained for each energy region by 
dividing the number of GG events in each bin by the to¬ 
tal number of CG events in that energy regions. At lower 
neutrino energies most of the GG events are selected in 
the downstream vertex-Z bin for the DS-escaping topol¬ 
ogy whereas at higher energies the DS-escaping events 
are distributed uniformly in vertex-Z. NonDS-escaping 
GG events are selected in all vertex-Z bins for low en¬ 
ergy neutrinos but higher energy neutrinos tend to be lo¬ 
cated in upstream vertex-Z bins. In addition, more high- 
energy neutrino events are selected in modules closer to 
the beam-axis. 


V. PROPAGATION OF SYSTEMATIC 
UNCERTAINTIES 


As described in Sec. IV F, the has terms with covari¬ 
ance matrices for systematic parameters. In this section. 
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FIG. 13. Probability density function (PDF) for each energy region defined in Sec. IV G: F^iy=0.5-0.8 GeV, 0.8-1.4 GeV, 
1.4-2.6 GeV, 2.6-4.0 GeV, >4.0 GeV. 
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we describe how the covariance matrices for the neutrino 
flux, the neutrino interaction model and the detector re¬ 
sponse, which were introduced in Secs. Ill A, IIIB and 
IV E, respectively, are constructed. 


A. Neutrino flux uncertainties 


The covariance matrix for each source of error on the 
neutrino flux perdiction, such as the horn current uncer¬ 
tainty, is calculated by taking the variation of the flux due 
to that error. The total covariance matrix is obtained by 
summing all these matrices and Figure 14 shows the cor¬ 
relation matrix derived from it. The energy binning in 
the covariance matrix is the same as that used to define 
the “global bin” in Table V. One can see that it is largely 
positively correlated. This correlation comes mainly from 
the uncertainty associated with hadron production at the 
target, which varies the neutrino flux in the same way at 
all INGRID modules. 


7 

6 

5 

4 

3 

2 

1 



37-42 

31-36 

25-30 

19-24 

13-18 

7-12 

1-6 


Energy 

bin 


Module 

group 


1-6 


1 


7-12 13-18 19-24 25-30 31-36 37-42 
3 4 5 6 


1 

0.8 

0.6 

0.4 


U 


where one of the NC systematic parameters has been 
changed by I a. The number of events is altered not only 
by the change in the cross section but also by changes 
in the event detection efficiency. The same normaliza¬ 
tion parameter, is used for all module groups, so the 
number of predicted events, becomes: 

^nc , £nc . ^nc (one ^ onc\ QiN 

39 h 39 ^>^3 ^39^' ^ ^ 

is estimated by combining NC events from all mod¬ 
ule groups. The fractional covariance for the topology 
bins is then calculated by varying each NC systematic 
parameter by Tier. 


=- 
^13 2 


_ ^nc, + lcr j^nc _ ^nc,+lcr ^ 




3 


JSfnc _ jync,-la jync _ jy_ 


nc, — la '' 






( 12 ) 


where is obtained by summing over all module 

groups: 
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We found that the total NC error is fully correlated be¬ 
tween all of the “NonDS-escaping” bins. Therefore, the 
7 “NonDS-escaping” bins are merged into a single bin. 
Figure 15 shows the correlation matrix for the NC inter¬ 
actions uncertainty, demonstrating that the event topol¬ 
ogy bins are almost fully correlated with each other. This 
correlation comes mainly from the NC normalization er¬ 
ror (see Table II). 


FIG. 14. Correlation matrix between module groups for the 
flux error. The energy binning used in this matrix is same as 
the “lobal bin” defined in Table V. 


B. Neutrino interaction model uncertainties 

Any systematic error in the CC interaction model 
would, by definition, alter the CC inclusive cross section 
itself. This is not the case for NC interaction model un¬ 
certainties, and so the systematic errors associated with 
these two processes are evaluated separately. 



VertexZ 1 2 3 4 5 6 7 

-escaping DS-escaping 

FIG. 15. Correlation matrix for the errors coming from 
NC interaction uncertainties. The binning on the y-axis is 
identical to that on the x-axis. 


7. Systematics uncertainty on NC interactions 

The uncertainty in the number of NC interactions in 
each bin is expressed by the normalization parameter: 

j^fne 

one ^ 39 

•^39 ~ IVfnc ’ 

39 

where is the predicted number of NC events for the 
topology and the module group. is the pre¬ 

dicted number of events in the same bin but for the case 


The total normalization error on the number of NC 
events is 27-30%, which is dominated by the NC nor¬ 
malization error, which is shown in Fig. 16. This gives a 
maximum error size of 5% on the total (CC+NC) number 
of events. 

2. Systematic uncertainty on CC interactions 

Varying the CC interaction parameters results in a 
change in both the neutrino cross section and the de- 
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FIG. 16. Fractional error from NC interaction uncertainties 
on the number of NC events. 


tection efficiency. In this analysis, for the systematic un¬ 
certainty on CC interactions, only the latter change is 
taken into account, since the result of the analysis will 
be the cross section itself. A Icr variation is applied to 
a CC interaction parameter and the new selection effi¬ 
ciency, is calculated. The change in the detection 
efficiency is then given by the ratio of the new efficiency 
to the nominal one: 


^cc ’ 


(14) 


where the indices i and j denote the energy bin and 
the topology bin, respectively. The predicted number 
of CC events is then modified using Wij\ 


i 

A fractional covariance between topology bins is then cal¬ 
culated for each CC interaction parameter using the mod¬ 
ified and nominal number of CC events. This is 

performed in the same way as for the NC interaction un¬ 
certainty, described earlier. The total covariance matrix 
is computed by summing up the individual matrix from 
each CC interaction parameter. The obtained correlation 
matrix and fractional error between the event topology 
bins are shown in Figs. 17 and 18. 


3. Systematic uncertainty on FSI 

For the pion FSI parameters, uncertainties on the ab¬ 
sorption, charged exchange, quasi-elastic, and inelastic 
scattering cross sections of the pion are taken into ac¬ 
count. These systematic errors are treated in a different 
to the previous interaction uncertainties because there 
are correlations between them. The INGRID data is fit¬ 
ted with obtained by changing each FSI parameter 
by Icr, and the difference between the fitted and the nom¬ 
inal result is taken as the systematic error. The effect of 
each FSI parameter on the measured cross section was 
found to be negligible, except for the pion absorption 



FIG. 17. Gorrelation matrix for the errors coming from 
GG interaction uncertainties. The binning on the y-axis is 
identical to that on the x-axis. 



NonDS-escaping DS-escaping 
FIG. 18. Fractional error from GG interaction uncertainties. 


uncertainty, which was then added in quadrature to the 
final result. 


C. Detector response uncertainties 


As described in Sec.IV E, each error source is catego¬ 
rized according to whether it produces correlation among 
topology bins or not. For the uncorrelated error sources, 
the iron mass and pileup correction, the individual sys¬ 
tematic errors are summed quadratically, and the total 
inserted into the denominator of the statistical term, 
{(j^det^ in Eq. (5)). 

3 9 

Eor the correlated error sources, the size of the error 
does not to vary between each module, so a covariance 
matrix is constructed from the topology bins using the 
average change over all modules. Namely, the errors are 
assumed to be fully correlated between module groups. 
Eor uncertainties from the event selection, the systematic 
error is evaluated by varying each selection threshold by 
Icr, and the resultant fractional variation in the number 
of selected events (= AN/N) computed for data and MC. 
The difference in AN/N between data and MC is then 
taken as the systematic error, calculated as: 


f ^Nobs \ _ f ANe^p \ 

V Nobs Jj V Ne.P Jj ’ 


( 16 ) 
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where the index j denotes the topology bin. If both a 
+lcr and — Icr variation are applied to the event selection 
then the following covariance is calculated: 

~ 2 + ^ 

If only a +1(7 change can be applied the the covariance 
becomes: 


y,, = A, • A, . (18) 

The statistical error of A^ is also calculated and added 
to Eq. (17) (or Eq. (18)). Einally, the total covariance is 
calculated by summing each individual covariance. Eig- 
ures 19 and 20 show the correlation matrix obtained and 
the size of the fractional error for each event topology 
bin. 



FIG. 19. Correlation matrix from the uncertainties in the 
detector response. The binning on the y-axis is identical to 
that on the x-axis. 



NonDS-escaping DS-escaping 

FIG. 20. Fractional error for the topology bins coming from 
uncertainties in the detector response. 


TABLE VI. Systematic error size on topology bins due to 
uncertainties in the detector response. 


Error type 

Error size (at maximum) 

Gorrelated error 

3% 

Uncorrelated error 

2% 


Table VI summarizes the size of the detector system¬ 
atic error for each error type. 


D. Uncertainty in pion multiplicities and 
secondary interactions 

Uncertainties associated with pion multiplicities and 
pion secondary interactions (SI) are treated in a differ¬ 
ent way to the systematics described above. These uncer¬ 
tainties are evaluated by comparing the underlying model 
with external data. Any observed difference is used to 
correct the nominal MC sample, and then the fit is 
performed using this “corrected” MC. The differences in 
the fitted values between the nominal and corrected MC 
are then taken as the systematic error on the final result. 

Pion multiplicity: 

In this analysis, the number of events is deter¬ 
mined from the number of reconstructed vertices, 
which are sometimes missed due to the pile-up of 
tracks from multiple neutrino interactions. Events 
with large numbers of tracks usually contain pi¬ 
ous, therefore the uncertainty associated with the 
pion multiplicity in these events needs to be con¬ 
sidered. This uncertainty is estimated by following 
the method described in Ref. [38]. The probabil¬ 
ity of an event having a pion multiplicity of n is 
expressed as: 


1 

P{z-, A, B,c) = — , (19) 

{n) = A + B\ogW^ , (20) 

where z=n/{n); {n) is the mean pion multiplicity 
and can be expressed by an approximate formula 
as in Eq. (20). The W is the hadronic invariant 
mass, and is expressed as: 

W'^ = {E^ + En - E^f - + pjv - Puf ( 21 ) 

with Eu^{p^^) and E^{p^) denoting the energy (3- 
momentum) of the and /i respectively; de¬ 
noting the Eermi momentum of the nucleon and 
the nucleon energy. and c are derived by fit¬ 

ting two external data sets [39, 40] with Eq. (19). 
These fitted parameters are compared with those 
used in NEUT, and the differences assigned as the 
systematic uncertainty. The parameters are used 
to produce a corrected MC sample which is input 
to the fif- The differences in the fitted values 
coming from the corrected and the nominal MC 
are taken as the systematic error on the final result 
due to pion multiplicity uncertainties. 

Pion SI: 

Hadrons produced in neutrino interactions can also 
interact whilst traveling through the detector, a 
process known as “secondary interaction (SI)”. In 
the INGRID simulation the pion SI processes are 
controlled by GEANT4. In order to evaluate the 
uncertainty in the pion SI model, the following in¬ 
teraction modes are considered: 
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TABLE VIL Summary of pion-nucleus scattering data used to evaluate the pion SI uncertainty. The reactive cross section is 
defined as the sum of all the inelastic cross sections. 



Hadrons Targets 

piab (MeV/c) 

Interaction type 

K. Nakai et al. [33] 

'K^ l'K~ 

Al, Ti, Cu, Sn, An 

83-395 

ABS 

D. Ashery et al. [34] 

l'K~ 

Li, C, 0, Al, Fe, Nb, Bi 

175-432 

Reactive, Elastic, QEL, ABS, SCX 

M.K. Jones et al. [35] 

TT^ 

C, Ni, Zr, Sn, Pb 

363-624 

QEL, ABS, SCX 

G.J. Gelderloos et al. [36] 

7r“ 

Li, C, Al, S, Ca, Cu, Zr, Sn, Pb 479-616 

Reactive 

B.W. Allardyce et al. [37] 

l'K~ 

C, Al, Ca, Ni, Sn, Ho, Pb 

710-2000 

Reactive 


Quasi-elastic scattering (QEL): 

The final state pion is the same type as the 
incoming pion. 

Absorption (ABS): 

The incident pion is absorbed by the nucleus, 
resulting in there being no pions in the final 
state. 

Single charge exchange (SCX): 

A charged pion interacts such that there is 
only one tt^ and no other pions in the final 
state. 


The existing experimental data used to evaluate 
this uncertainty is summarized in Table VII. In 
the table, the reactive cross section is defined as 
c^totai - CTeiastic, whcrc (Jtotai IS the total cross sec¬ 
tion and cTeiastic is the elastic cross section. As seen 
in the table, D. Ashery et al provides various cross 
sections across a range of pion momenta and tar¬ 
get nuclei, including iron. The other data do not 
include measurements on iron. For these data, an 
A-dependent scaling is applied in order to extract 
the cross section on iron. To evaluate the system¬ 
atic uncertainty coming from pion SI, we first tune 
the pion cross section in the momentum region cov¬ 
ered by the data in Table VII. Second, for the lower 
energy region not covered by data, the ABS 7r+(7r“) 
cross section <20(30) MeV is kept constant, moti¬ 
vated by the microscopic calculation from [41]. The 
QEL cross section is extrapolated to 0 at 0 MeV. 
For the higher energy region, each of the cross sec¬ 
tions above is tuned on the basis that the size of the 
reactive cross section is conserved, since the cross 
sections predicted by GEANT4 are in agreement 
with the experimental data in this energy region. 
This study gives 4 corrected MC samples in total, 
each of which is then used to fit the T2K data. 
Finally, the size of the systematic error on the fk 
parameters due to the uncertainty on pion FSI is 
calculated as follows: 


\ i=i 


( 22 ) 


where the index i denotes the corrected MC 


sample and fk is the fitted normalization parameter 
for the energy bin. 


VI. RESULT 

In this section, we present the result of this inclusive 
CC cross-section measurement. Section VIA shows the 
data set used in this analysis, whilst Secs. VIB and VIC 
describe the output from the give a summary 

of this result. 


A. Data set 

Figure 21 shows the observed and predicted topology 
distributions in all module groups for the data set used 
in this analysis, corresponding to 6.27 x 10^^ p.o.t. The 
number of observed events for the NonDS-escaping topol¬ 
ogy is 3-10% smaller than expected. 


B. Cross section fit 

The topology distributions after the data fit are shown 
in Fig. 22. As seen in the figure, the predicted topology 
distributions with the fitted cross-section normalization 
parameters applied agree well with the observed distri¬ 
butions. 

Table VIII shows the cross-section normalization pa¬ 
rameters, fi = l-\- Afi (i=0-4), obtained by fitting the 
INGRID data, where the A/^’s are those from Eq. (6). 

The fitted values for the flux, detector, CC interaction, 
and NC interaction uncertainty parameters are shown 
in Fig. 23. A large deviation from 0 is seen for all the 
NC interaction uncertainty parameters. As described in 
Sec. VIA, the prediction overestimates the number of 
NonDS-escaping events by 3-10%. 

The fitter preferentially reduces the number of NC 
events to match the predicted topology distribution to 
the observed one. Since the NC interaction uncertainty 
parameters are almost fully correlated among topologies, 
as shown in Fig. 15, all the parameters move toward neg¬ 
ative values. There are jumps seen in the CC interaction 
and detector uncertainty parameters. Both of the jumps 
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FIG. 21. Topology distribution for NonDS-escaping (top) and DS-escaping (bottom) for all module groups. The predicted 
events, before the ht, are categorized as CC events, NC events, and either intrinsic beam + Ue backgrounds or backgrounds 
from the wall, and shown as a stacked histogram. The CC events are subdivided into 6 true neutrino energy regions. 
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FIG. 22. Topology distribution for NonDS-escaping (top) and DS-escaping (bottom) events for all module groups after the 
data fit. The predicted events are categorized as CC events, NC events, and either intrinsic beam backgrounds or 

backgrounds from the wall, and shown as a stacked histogram. The CC events are subdivided into 6 true neutrino energy 
regions. 
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TABLE VIIL Fitted values of the cross-section normalization 
parameters, /i = 1 + Afi. 


Fit parameter 

Fit result 

/o (0.5 GeV) 

3.560 ± 0.508 

/i (0.8 GeV) 

0.637 ± 0.180 

/2 (1.4 GeV) 

1.324 ± 0.181 

h (2.6 GeV) 

0.800 ± 0.211 

U (4.0 GeV) 

0.712 ± 0.120 


155.4 


cross-section normalization at 2.0 GeV is less affected by 
most of the systematic uncertainties than the other nor¬ 
malizations, as shown in Table X. The reason for this is 
as follows. In Fig. 13, one can see that the PDF is well 
differentiated for both module group and event topology 
at E^=lA-2.6 GeV (/ 2 -/ 3 ) and E^=2.6A.O GeV {fs-h), 
while poor differentiation is seen for E^ >4.0 GeV (/ 4 ), 
which results in weak correlation with the other cross- 
section normalizations. Thus, sensitivity to the cross- 
section normalization is good for Ey=2A GeV (= 
but limited for Ey=3.3 GeV (= 
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FIG. 23. Fitted values of the cross-section normalization 
(top) and systematic (bottom) parameters. 


appear at the boundary between parameters for NonDS- 
escaping and DS-escaping events. 

In order to derive the normalization factor for the cross 
section, we take the average of the neighboring fitted fi 
parameters. The obtained cross-section normalizations 
are: 

/(l.l GeV) = 

= 0.980 ±0.115 , 

/(2.0 GeV) = 'krA 

= 1.062 ±0.123 , 

/(3.3 GeV) = 

= 0.756 ±0.136 , 


TABLE IX. Contribution to the uncertainty on the htted pa¬ 
rameters (/0-/4) from each error source. 


Error source 

fo 


h 

fs 

h 


(0.5 GeV) (0.8 GeV) (1.4 GeV) (2.6 GeV) (4.0 GeV) 

Statistical error 

18.7% 

6 . 0 % 

4.5% 

4.8% 

1.4% 

Flux + Stat. 

26.0% 

7 . 9 % 

12.8% 

14.5% 

9.3% 

Detector + Stat. 

33.8% 

10 . 0 % 

7.2% 

7.0% 

3.0% 

Interaction (cc) + Stat. 

30.6% 

9 . 3 % 

6.8% 

7.2% 

3.8% 

Interaction (nc) + Stat. 

22.6% 

6.6% 

6.4% 

5.8% 

2.0% 

pion ESI 

±0.4% 

-0.2% 

±0.3% 

-0.8% 

±2.0% 

-3.3% 

±4.0% 

-3.6% 

±3.5% 

-5.2% 

pion multiplicity 

2.6% 

8.8% 

0.7% 

12.4% 

9.4% 

pion SI 

7.3% 

9.5% 

9.4% 

11.4% 

2.9% 


TABLE X. Contribution to the uncertainty on the cross- 
section normalization at 1.1, 2.0, and 3.0 GeV from each error 
source. 


Error source 

1.1 GeV 

2.0 GeV 

3.3 GeV 

Statistical error 

2 . 0 % 

0 . 6 % 

2 . 3 % 

Flux + Stat. 

7 . 6 % 

9 . 0 % 

8 . 4 % 

Detector + Stat. 

4 . 3 % 

0 . 9 % 

3 . 9 % 

Interaction (cc) + Stat. 

3 . 7 % 

0 . 8 % 

4 . 8 % 

Interaction (nc) + Stat. 

2 . 4 % 

0 . 9 % 

3 . 2 % 

pion FSI 

+1.0% 

-1.9% 

0 . 5 % 

+3.7% 

-2.9% 

pion multiplicity 

3 . 3 % 

5 . 1 % 

2 . 1 % 

pion SI 

5 . 6 % 

2 . 0 % 

6 . 9 % 


Figures 24 and 25 show the error and correlation matri¬ 
ces for the 5 fitted parameters (A/ 0 -A/ 4 ) and the cross- 
section normalizations at 1 . 1 , 2 . 0 , and 3.3 GeV, respec¬ 
tively. 


Table IX summarizes the uncertainty on the fitted 
cross-section normalization parameters, broken down by 
error source. The errors on the combined normalization 
parameters are summarized in Table X. The largest sys¬ 
tematic error source is the flux uncertainty, which gives 
a 8-9% error on the cross-section normalization. The 


C. Summary 

In the previous section, five individual fitting parame¬ 
ters were extracted with the least method, and used 
to calculate the following measured cross sections at en- 
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FIG. 24. Error (left) and correlation (right) matrices for the 
5 fitted parameters. In both of the matrices, the binning on 
the y-axis is identical to that on the x-axis. 


and agrees with the MINOS data point and GENIE at 
Ej^=3.3 GeV. Further reduction of the systematic error 
size would help in differentiating the neutrino models at 
the higher energy transition. This reduction could be 
made by utilizing neutrino beams covering a wider range 
of off-axis angles, which would provide a “model inde¬ 
pendent” way to infer the neutrino energy. This analy¬ 
sis demonstrates the feasibility of using different off-axis 
samples from the same neutrino beam to measure the 
energy dependence of neutrino interactions, which will 
provide useful information for future neutrino oscillation 
analyses. 



FIG. 25. Error (left) and correlation (right) matrices for the 
cross-section normalization at 1.1, 2.0, and 3.3 GeV. In both 
of the matrices, the binning on the y-axis is identical to that 
on the x-axis. 


ergies of 1.1, 2.0, and 3.3 GeV: 

cr^^(l.l GeV) = 1.10 ±0.15 (10“^^cm^/nucleon) , 
cr^^(2.0 GeV) = 2.07 ±0.27 (10“^^cmVnucleon) , 
cr^^(3.3 GeV) = 2.29 ± 0.45 (10“^^cmVnucleon) . 


VII. CONCLUSION 

In this paper, we have reported the measurement of 
the energy dependent inclusive charged current cross 
section on iron using the T2K INGRID detector and the 
T2K neutrino beam. The unique variation in the neu¬ 
trino flux across the INGRID modules, along with event 
kinematic information, was used to produce event sam¬ 
ples sensitive to neutrinos with energies from 1-3 GeV. 
These were used to extract the inclusive CC muon neu¬ 
trino interaction cross section on iron at 1.1, 2.0 and 3.3 
GeV, using data corresponding to 6.27 x 10^^ p.o.t. This 
result is consistent with the predictions of the NEUT and 
GENIE neutrino interaction generators. This is the first 
measurement to use the off-axis effect to measure neu¬ 
trino cross sections as a function of energy and demon¬ 
strates the potential of this technique to provide useful 
information for future neutrino experiments. 


Figure 26 shows these results compared to other measure¬ 
ments [1, 2] and the neutrino event generators, NEUT 
(ver.5.1.4.2) and GENIE (ver.2.8.0). These measure¬ 
ments are consistent with the energy dependent cross 
section measured by the MINOS near detector and the 
previous, flux-averaged, cross section measured by IN¬ 
GRID, which used a subset of the data included in this 
analysis. 

This analysis utilizes the different off-axis angle tech¬ 
nique together with the final state kinematics of the out¬ 
going lepton to enhance sensitivity to the neutrino en¬ 
ergy. Using final state lepton kinematics in this analysis 
makes the result sensitive to uncertainties in the neu¬ 
trino interaction model, which increases the uncertainty 
on the final measurement. As a result, the errors achieved 
are not small enough to distinguish between the neutrino 
models used in the different event generators. Neverthe¬ 
less, the result seems to prefer NEUT for Ejy <2 GeV 
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